Spin-orbital liquids in non-Kramers magnet on Kagome lattice 
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Localized magnetic moments with crystal-field doublet or pseudo-spin 1/2 may arise in correlated insulators 
with even number of electrons and strong spin-orbit coupling. Such a non-Kramers pseudo-spin 1/2 is the 
consequence of crystalline symmetries as opposed to the Kramers doublet arising from time-reversal invariance, 
and is necessarily a composite of spin and orbital degrees of freedom. We investigate possible spin-orbital 
liquids with fermionic spinons for such non-Kramers pseudo-spin 1/2 systems on the Kagome lattice. Using the 
projective symmetry group analysis, we find ten new phases that are not allowed in the corresponding Kramers 
systems. These new phases are allowed due to unusual action of the time reversal operation on non-Kramers 
pseudo-spins. We compute the spin-spin dynamic structure factor that shows characteristic features of these 
non-Kramers spin-orbital liquids arising from their unusual coupling to neutrons, which is therefore relevant for 
neutron scattering experiments. We also point out possible anomalous broadening of Raman scattering intensity 
that may serve as a signature experimental feature for gapless non-Kramers spin-orbital liquids. 
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I. INTRODUCTION 

The low energy magnetic degrees of freedom of a Mott 
insulator, in the presence of strong spin-orbit coupling, are 
described by states with entangled spin and orbital wave 
functions.''^ In certain crystalline materials, for ions with 
even numbers of electrons, a low energy spin-orbit entan- 
gled "pseudo-spin"- 1/2 may emerge, which is not protected 
by time-reversal symmetry (Kramers degeneracy)'^ but rather 
by the crystal symmetries.-'"* Various phases of such non- 
Kramers pseudo-spin systems on geometrically frustrated lat- 
tices, particularly various quantum paramagnetic phases, are 
of much recent theoretical and experimental interest in the 
context of a number of rare earth materials including frus- 
trated pyrochlores^"^ and heavy fermion systems."*-" 

In this paper, we explore novel spin-orbital liquids that may 
emerge in these systems due to the unusual transformation of 
the non-Kramers pseudo-spins under the time reversal trans- 
formation. Contrary to Kramers spin- 1/2, where the spins 
transform as S — > — S under time reversal,'^ here only one 
component of the pseudo-spin operators changes sign under 
time reversal: {a^,a^,a^} — > {ct^ct^, —tr^}.-'''* This is be- 
cause, due to the nature of the wave-function content, the 
(73 component of the pseudo-spin carries a dipolar magnetic 
moment while the other two components carry quadrupolar 
moments of the underlying electrons. Hence the time rever- 
sal operator for the non-Kramers pseudo-spins is given by 
T = a^K (where K is the complex conjugation operator), 
which allows for new spin-orbital liquid phases. Since the 
magnetic degrees of freedom are composed out of wave func- 
tions with entangled spin and orbital components, we prefer 
to refer the above quantum paramagnetic states as spin-orbital 
liquids, rather than spin liquids. 

Since the degeneracy of the non-Kramers doublet is pro- 
tected by crystal symmetries, the transformation properties of 
the pseudo-spin under various lattice symmetries intimately 



depend on the content of the wave-functions that make up 
the doublet. To this end, we focus our attention on the ex- 
ample of Praseodymium ions (Pr''+) in a local D^j^ environ- 
ment, which is a well known non-Kramers ion that occurs in 
a number of materials with interesting properties.''"^ Such an 
environment typically occurs in Praseodymium pyrochlores 
given by the generic formulae Pr2TM207, where TM(= Zr, 
Sn, Hf, or Ir) is a transition metal. In these compounds, the 
Pr-^^ ions host a pair of 4/ electrons which form a J = 4 
ground state manifold with S — 1 and i = 5, as expected due 
to Hund's rules. In terms of this local environment we have 
a nine fold degeneracy of the electronic states. "* This degener- 
acy is broken by the crystalline electric field. The oxygen and 
TM ions form a D^d local symmetry environment around the 
Pr-^^ ions, splitting the nine fold degeneracy. A standard anal- 
ysis of the symmetries of this system (see appendix A) shows 
that the J = 4 manifold splits into three doublets and three 



singlets (Tj 



3Eg + 2Aig + A2g) out of which one of the 



doublets is found to have the lowest energy, usually well sepa- 
rated from the other crystal field states.^ This doublet (details 
in Appendix A), formed out of a linear combination of the 
J^ ~ ±4 with ,Jz = ±1 and J^ = ±2 states, is given by 



l±) 



aim 



±4) ±/3|to = ±1) -7|to = =F2). (1) 



The non-Kramers nature of this doublet is evident from the 
nature of the "spin" raising and loweiing operators within 
the doublet manifold; the projection of the angular momen- 
tum raising and loweiing operators to the space of doublets 
is zero (VJ^Vlcr) = where V projects into the doublet 
manifold). However, the projection of the J^ operator to this 
manifold is non-zero, and describes the z component of the 
pseudo-spin (cr^). In addition, there is a non-trivial projection 
of the quadrupole operators { J^ , J^ } in this manifold. These 
have off-diagonal matrix elements, and are identified with the 
pseudo-spin raising and loweiing operators (cr* = a^ ± ia"^). 
In a pyrochlore lattice the local D^d axes point to the centre 
of the tetrahedra.* On looking at the pyrochlore lattice along 



the [111] direction, it is found to be made out of alternate 
layers of Kagome and triangular lattices. For each Kagome 
layer (shown in Fig. 1) the local Ds^ axes make an angle of 
cos~'^{^y2/3) with the plane of the Kagome layer. We imag- 
ine replacing the Pr'^+ ions from the triangular lattice layer 
with non-magnetic ions so as to obtain decoupled Kagome 
layers with Pr'^^ ions on the sites. The resulting structure is 
obtained in the same spirit as the now well-known Kagome 
compound Herbertsmithite was envisioned. As long as the lo- 
cal crystal field has D^d symmetry, the doublet remains well 
defined. A suitable candidate non-magnetic ion may be iso- 
valent but non-magnetic La'^^. Notice that the most extended 
orbitals in both cases are the fifth shell orbitals and the crystal 
field at each Pr'^^" site is mainly determined by the surrounding 
oxygens and the transition metal element. Hence, we expect 
that the splitting of the non-Kramers doublet due to the above 
substitution would be very small and the doublet will remain 
well defined. In this work we shall consider such a Kagome 
lattice layer and analyze possible Z2 spin-orbital liquids, with 
gapped or gapless fermionic spinons. 

The rest of the paper is organized as follows. In Sec. 11, we 
begin with a discussion of the symmetries of the non-Kramers 
system on a Kagome lattice and write down the most general 
pseudo-spin model with pseudo-spin exchange interactions up 
to second nearest neighbours. In Sec. Ill formulate the pro- 
jective symmetry group (PSG) analysis for singlet and triplet 
decouplings. Using this we demonstrate that the non-Kramers 
transformation of our pseudo-spin degrees of freedom under 
time reversal leads to a set of ten spin-orbital liquids which 
cannot be realized in the Kramers case. In Sec. IV we de- 
rive the dynamic spin-spin structure factor for a representa- 
tive spin liquid for the case of both Kramers and non-Kramers 
doublets, demonstrating that experimentally measurable prop- 
erties of these two types of spin-orbital liquids differ qualita- 
tively. Finally, in Sec. V, we discuss our results, and propose 
an experimental test which can detect a non-Kramers spin- 
orbital liquid. The details of various calculations are discussed 
in different appendices. 



11. SYMMETRIES AND THE PSEUDO-SPIN 
HAMILTONIAN 

Since the local D^d axes of the three sites in the Kagome 
unit cell differ from each other a general pseudo-spin Hamil- 
tonian is not symmetric under continuous global pseudo-spin 
rotations. However, it is symmetric under various symmetry 
transformations of the Kagome lattice as well as time rever- 
sal symmetry. Such symmetry transformations play a major 
role in the remainder of our analysis. We start by describing 
the effect of various lattice symmetry transformations on the 
non-Kramers doublet. 

We consider the symmetry operations that generate the 
space group of the above Kagome lattice. These are (as shown 
in Fig. 2(a)) 

• Ti,T2 : generate the two lattice translations. 

• a — C2I '■ (not to be confused with the pseudo-spin 




FIG. 1 : A Kagome layer, in the pyrochlore lattice environment. We 
consider sites labelled z and z' replaced by non-magnetic ions, de- 
coupling the Kagome layers. The local axis at the u,v and w sites 
point towards the center of the tetrahedron on which these lie. 




FIG. 2: (color online) (a) The symmetries of the Kagome lattice. 
Also shown are the labels for the sublattices and the orientation of the 
local z-axis. (b) Nearest and next nearest neighbour bonds. Colors 
refer to the phases (f)r^r' and (j>'^ ^, , with these being on blue bonds, 
1 on green bonds and 2 on red bonds. 



operators which come with a superscript) where / is the 
three dimensional inversion operator about a plaquette 
center and C2 refers to a two-fold rotation about a line 
joining two opposite sites on the plaquette. 

• Sq ^ C^I : where C3 is the threefold rotation op- 
erator about the center of a hexagonal plaquette of the 
Kagome lattice. 

• T = (JiK : Time reversal. 

Here, we consider a three dimensional inversion operator 
since the local D^d axes point out of the Kagome plane. The 
above symmetries act non-trivially on the pseudo-spin degrees 
of freedom, as well as the lattice degrees of freedom. The ac- 
tion of the symmetry transformations on the pseudo-spin op- 
erators is given by. 



Se:{<y\<J+,<j-}^{a\ 



Loa , u)u 



-}, 



T : {c7^,a+,a-} -^ {-cr^, ff-, cr+}, 
C^:{a\a+,a-}^{-a\<j-,a+}, 
Ti:{a^,a+,a-}^{a^,a+,a-}, 
T2 : {cr\a+,<j-} -^ {a\a+,a-}, 



(2) 



(uj = Lu^ = e*T"). Operationally their action on the doublet 
(1+) K)) '^an be written in form of 2 x 2 matrices. The 
translations Ti , T2 act trivially on the pseudo-spin degrees of 



freedom, and the remaining operators act as 

T^aiK, cr = o-i, 56 = 



H) 

uj 



(3) 



where K refers to complex conjugation. The above expres- 
sions can be derived by examining the effect of these operators 
on the wave-function describing the doublet (Eq. 1). 

We can now write down the most generic pseudo-spin 
Hamiltonian allowed by the above lattice symmetries that is 
bilinear in pseudo-spin operators. The form of the time- 
reversal symmetry restricts our attention to those products 
which are formed by a pair of a^ operators or those which 
mix the pseudo-spin raising and lowering operators. Any 
term which mixes a'^ and cr* changes sign under the sym- 
metry, and can thus be excluded. Under the C3 transforma- 
tion about a site, the terms C3 : a'^af, — > a^ (r)'^c (r>) ^^'^ 
C3 ■ cr+a~, -^ a'^^^^^a~_^^^,y However, the term cr+a+ (and 
its Hermitian conjugate) gain additional phase factors when 
transformed; under the C3 symmetry transformation, this term 
becomes C3 : ct+ct^ — > tJCTp (r)'^c (r')- ^" addition, un- 
der the (7 symmetry, this term transforms as a : a^a^, ^■ 
cr^i^)0'^/j.,\- Thus the Hamiltonian with spin-spin exchange 
interactions up to next-nearest neighbour is given by 



Heff = Jnn ^ [^r^Cr^, + 2{6a+a^, + h.C.) 

{r.r') 

+2q{e r^cr+cr+ + h.c.)] 

+ Jnnn ^ [ct^ct^/ + 2{S' a^^ a^^, + h.c.) 

((r,r')> 

+2q'{e'-^^a+a++h.c.)], (4) 

where 4> and (p' take values 0, 1 and 2 depending on the bonds 
on which they are defined (Fig. 2(b)). 



III. SPINON REPRESENTATION OF THE PSEUDO-SPINS 
AND PSG ANALYSIS 

Having written down the pseudo-spin Hamiltonian, we now 
discuss the possible spin-orbital liquid phases. We do this in 
two stages in the following sub-sections. 



A. Slave fermion representation and spinon decoupling 

In order to understand these phases, we will use the 
fermionic slave-particle decomposition of the pseudo-spin op- 
erators. At this point, we note that the pseudo-spins satisfy 
S = 1/2 representations of a "SU(2)" algebra among their 
generators (not to be confused with the regular spin rotation 
symmetry). We represent the pseudo-spin degrees of freedom 
in terms of a fermion bilinear. This is very similar to usual 
slave fermion construction for spin liquids'^''*. We take 



2/la[^^]"/3/j73' 



(5) 



where a, /? =t, i is defined along the local z axis and /^ (/) 
is an 5 = 1/2 fermionic creation (annihilation) operator. Fol- 
lowing standard nomenclature, we refer to the /(/^) as the 
spinon annihilation (creation) operator, and note that these 
satisfy standard fermionic anti-commutation relations. The 
above spinon representation, along with the single occupancy 
constraint 



fl^M + flu = 1, 



(6) 



form a faithful representation of the pseudo-spin- 1/2 Hilbert 
space. The above representation of the pseudo-spins, when 
used in Eq. 4, leads to a quartic spinon Hamiltonian. Fol- 
lowing standard procedure,'^^''* this is then decomposed us- 
ing auxiliary fields into a quadratic spinon Hamiltonian (af- 
ter writing down the corresponding Eucledian action). The 
mean field description of the phases is then characterized by 
the possible saddle point values of the auxiliary fields. There 
are eight such auxiliary fields per bond, corresponding to 



Xtj = ifLfja)*; mj = {f^a [ir^] fjp) 



(7a) 
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-'-1.,/'/^) 



(7b) 



where t" (a = 1, 2, 3) are the Pauli matrices. While Eq. 7a 
represents the usual singlet spinon hopping (particle-hole) and 
pairing (particle-particle) channels, Eq. 7b represents the cor- 
responding triplet decoupling channels. Since the Hamilto- 
nian (Eq. 4) does not have pseudo-spin rotation symmetry, 
both the singlet and the triplet decouplings are necessary.""''^ 
From this decoupling, we obtain a mean-field Hamiltonian 
which is quadratic in the spinon operators. We write this com- 
pactly in the following form'^ (subject to the constraint Eq.6) 



Ho - 2^Ji]fi Uijfj, 



(8) 
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— fi-f fii fil fit 
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(9) 


Urj 


= e^fs"r^, 


(10) 


E" 


= p"0/; T^^I^T^, 


(11) 



where p" are the Identity (for a = 0) and Pauli matrices 
(a = 1, 2, 3) acting on pseudo-spin degrees of freedom, and 
r" represents the same in the gauge space. We immediately 
note that 



[E",r''] =0 Va,/3. 



(12) 



The requirement that our Hq be Hermitian restricts the coeffi- 
cients S,ij to satisfy 






(13) 



for a, 6 G {1,2,3}. The relations between f^jS and 
{Xij: Vij^ '^ij- ■'-^ij} ^^ given in Appendix C.'^ As a straight 
forward extension of the SU{2) gauge theory formulation for 



spin liquids, '^''^ we find that Hq is invariant under the gauge 
transformation 
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Z^lj"j^ 



(14) 
(15) 



where the Wi matrices are SU(2) matrices of the form Wi = 
gifo- (f = (r\r2,r3)). Noting that the physical pseudo- 
spin operators are given by 



1 ^- - 



(16) 



Eq. 12 shows that the spin operators, as expected, are gauge 
invariant. It is useful to define the "E-components" of the Uij 
matrices as follows: 



where 



and 
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(17) 



(18) 



(19) 



(20) 



Under global spin rotations the fermions transform as 

where V is an SU(2) matrix of the form V = e'^^ (E = 
{E^ , E^ , E'^}). So while V? (the singlet hopping and pairing) 
is invariant under spin rotation, {V/,, Vf,, Vf,} transforms as 
a vector as expected since they represent triplet hopping and 
pairing amplitudes. 



B. PSG Classification 

We now classify the non-Kramers spin-orbital liquids based 
on projective representation similar to that of the conven- 
tional quantum spin liquids.'^ Each spin-orbital liquid ground 
state of the quadratic Hamiltonian (Eq 11) is character- 
ized by the mean field parameters (eight on each bond, 
X, rj, E^,E'^,E^, D'^,D'^,D^, or equivalently Uij). However, 
due to the gauge redundancy of the spinon parametrization 
(as shown in Eq. 15), a general mean-field ansatz need not 
be invariant under the symmetry transformations on their own 
but may be transformed to a gauge equivalent form without 
breaking the symmetry. Therefore, we must consider its trans- 
formation properties under a projective representation of the 
symmetry group. '^^ For this, we need to know the various pro- 
jective representations of the lattice symmetries of the Hamil- 
tonian (Eq. 4) in order to classify different spin-orbital liquid 
states. 

Operationally, we need to find different possible sets of 
gauge transformations {GG} which act in combination with 



the symmetry transformations {SG} such that the mean-field 
ansatz Uij is invariant under such a combined transformation. 
In the case of spin rotation invariant spin-liquids (where only 
the singlet channels x and 77 are present), the above statement 
is equivalent to demanding the following invariance: 

[/,, = [GsS] U,, [GsS\^ = Gs{i)Us(i)SU)Gl{]). (21) 

where S E SG is a symmetry transformation and Gs G GG 
is the corresponding gauge transformation. The different pos- 
sible {Gs|V5 e SG} give the possible algebraic PSGs that 
can characterize the different spin-orbital liquid phases. To 
obtain the different PSGs, we start with various lattice sym- 
metries of the Hamiltonian. The action of various lattice 
transformations'^ is given by 



Ti 



ix,y,s) -^ {x + l,y,s); 



T2 ■■ix,y,s) -> {x,y + l,s); 
a ■.{x,y,u) -^ {y,x,u); 

(x,y,v) -^ {y,x,w); 

(x,y,w) -^ {y,x,v); 
Se -{x, y,u) ^ {-y-l,x + y + 1, v); 

{x,y,v) ^- {-y,x + y,w); 

{x, y, w) -> {~y-l,x + y, u); 



(22) 



where {x, y) denotes the lattice coordinates and s e {u, v, w} 
denotes the sub-lattice index (see figure 2). 

In terms of the symmetries of the Kagome lattice, these op- 
erators obey the following conditions 



T- = a- = (iej" = e, 



(^6^^ 



T-'gT = e V.9 G SG, 

a'^T^^aT2 = e, 
a^^T2^aTi = e, 

Sq^T2^TiS6T2 = e, 
a^^SeaSg = e. 



(23) 



In addition, these commutation relations are valid in terms of 
the operations on the pseudo-spin degrees of freedom, as can 
be verified from Eq. 3. 

In addition to the conditions in Eq. 23, the Hamiltonian is 
trivially invariant under the identity transformation. The in- 
variant gauge group (IGG) of an ansatz is defined as the set of 



all pure gauge transformations Gj such that Gj : Ui^ 



U,, 



The nature of such pure gauge transformations immediately 
dictates the nature of the low energy fluctuations about the 
mean field state. If these fluctuations do not destabilize the 
mean-field state, we get stable spin liquid phases whose low 
energy properties are controlled by the IGG. Accordingly, 
spin liquids obtained within projective classification are pri- 
marily labelled by their IGGs and we have Z2,U{1) and 
SU{2) spin liquids corresponding to IGGs of Z2, U{1) and 



SU{2) respectively. In this work we concentrate on the set of 
Z2 "spin Hquids" (spin-orbital liquids with a Z2 IGG). 

We now focus on the PSG classification. As shown in Eq. 
2, in the present case, the pseudo-spins transform non-trivially 
under different lattice symmetry transformations. Due to the 
presence of the triplet decoupling channels the non-Kramers 
doublet transforms non-trivially under lattice symmetries (Eq. 
3). Thus, the invariance condition on the UijS is not given by 
Eq. 21, but by a more general condition 



C/., = [GsS] [/,, [GsS]^ = Gsi^)<l>s [Usv)Sif)] gUj). 



Here 



(l^s [Us{i)s{j)] = T^sUs(i)s(3)T^. 



(24) 



(25) 



and Vs generates the pseudo-spin rotation associated with the 
symmetry transformation {S) on the doublet. The matrices 
Vs have the form 



Vs 



2 



Va =X't =i*E\ 



iV3. 



Vt, = Vt, = S°. 



(26) 
(27) 



Under these constraints, we must determine the relations 
between the gauge transformation matrices Gs{i) for our set 
of ansatz. The additional spin transformation (Eq. 25) does 
not affect the structure of the gauge transformations, as the 
gauge and spin portions of our ansatz are naturally separate 
(Eq. 12). In particular, we can choose to define our gauge 
transformations such that 

Gs : U.J = Gs : ^f E"r^ ^ C/S"G^(*)r'^G5(j), 

(28) 

S:U,,^S: 4f E"r^ ^ Csf,)so)^sS"2?Ir^ (29) 

where we have used the notation Gs '■ Utj = Gg{i)UijGs{j) 
and so forth. As a result, we can build on the general con- 
struction of Lu et al. ^^ to derive the form of the gauge trans- 
formation matrices. The details are given in Appendix B. 

A major difference arises when examining the set of alge- 
braic PSGs for Z2 spin liquids found on the Kagome lattice 
due to the difference between the structure of the time re- 
versal symmetry operation on the Kramers and non-Kramers 
pseudo-spin-l/2s. In the present case, we find there are 30 
invariant PSGs leading to thirty possible spin-orbital liquids. 
This is in contrast with the Kramers case analysed by Lu et 
al.,^^ where ten of the algebraic PSGs cannot be realized as 
invariant PSGs, as all bonds in these ansatz are predicted to 
vanish identically due to the form of the time reversal oper- 
ator, and hence there are only twenty possible spin liquids. 
However, with the inclusion of spin triplet terms and the non- 
Kramers form of our time reversal operator, these ansatz are 
now realizable as invariant PSGs as well. The time reversal 
operator, as defined in Appendix B, acts as 



T : ef/S° 



^ Ctf^'^T^, 



where f"'' = ^"'^ if a G {1, 2} and 1"^^ = -^"'^ if a G {0, 3}. 
The projective implementation of the time-reversal symmetry 
condition (Eq. 23) takes the form (see Appendix B) 



[GT{i)f^VTl Vz, 



(31) 



where Grii) is the gauge transformation associated with time 
reversal operation and 77^ = ±1 for a Z2 IGG. 

Therefore, the terms allowed by the time reversal symmetry 
to be non zero are, for tjt — 1, 



e',e\e\e\e\e\e\e\ 



(32) 



and for tjt = — 1, with the choice Grii) — iT^ (see appendix 
B), 



t02 ^03 flO <jll ^20 ^21 <j32 /;33 



(33) 



This contrasts with the case of Kramers doublets, in which no 
terms are allowed for 777- = —1, and for 777- = —1 the allowed 
terms are 



/:02 /:03 /:12 fl3 /:22 /:23 ^32 /:33 



(34) 



Further restrictions on the allowed terms on each link arise 
from the form of the gauge transformations defined for the 
symmetry transformations. All nearest neighbour bonds can 
then be generated from C/y defined on a single bond, by per- 
forming appropriate symmetry operations. 

Using the methods outlined in earlier works (Ref. 13, 15) 
we find the minimum set of parameters required to stabilize Z2 
spin-orbital liquids. We take into consideration up to second 
neighbour hopping and pairing amplitudes (both singlet and 
triplet channels). The results are listed in Table I. 

The spin-orbital liquids listed from 21 — 30 are not allowed 
in the case of Kramers doublets and, as pointed out before, 
their existence is solely due to the unusual action of the time- 
reversal symmetry operator on the non-Kramers spins. Hence 
these ten spin-orbital liquids are qualitatively new phases that 
may appear in these systems. Of these ten phases, only two 
(labelled as 21 and 22 in Table I) require next nearest neigh- 
bour amplitudes to obtain a Z2 spin-orbital liquid. For the 
other eight, nearest neighbour amplitudes are already suffi- 
cient to stabilize a Z2 spin-orbital liquid. 

It is interesting to note (see below) that bond-pseudo-spin- 
nematic order (Eq. 35 and Eq. 36) can signal spontaneous 
time-reversal symmetry breaking. Generally, since the triplet 
decouplings are present, the bond nematic order parameter for 
the pseudo-spins^''^^ 

Q,f = {[S^SI + 5f 5;) /2 - r^(S^, • Sj)li), (35) 

as well as vector chirality order 



Jij — \^i X Jj), 



(36) 



(30) 



are non zero. Since the underlying Hamiltonian Eq. 4) gener- 
ally does not have pseudo-spin rotation symmetry, the above 
non-zero expectation values do not spontaneously break any 



TABLE I: Symmetry allowed terms: We list the terms allowed to be non-zero by symmetry, for the 30 PSGs determined by Yuan-Ming Lu et 
a/'^. The PSGs listed together are those with 7712 = ±1 and all other factors equal. Included are terms allowed on nearest and next-nearest 
neighbour bonds, as well as chemical potential terms V which can be non zero on all sites for certain spin-orbital liquids. Also included is the 
distance of bond up to which we must include in order to gap out the gauge fluctuations to Z2 via the Anderson-Higgs mechanism '^^. Only 
PSGs 9 and 10 can not host Z2 spin-orbital liquids with up to second nearest neighbour bonds. 



No. 


A. 


n.n. 


n.n.n. 


Z2 


1-2 


p2 p3 


£lO £21 £02 £03 £32 £33 


£lO £21 £02 £03 £32 £33 


n.n. 


3-4 





flO c21 c02 c03 c32 c33 
S iS iS iS iS iS 


e\e' 


n.n. 


5-6 


p3 


£lO £21 £02 £03 £32 £33 


£lO £21 £03 £33 


n.n. 


7-8 





e\e' 


e\e\e\e\e\e' 


n.n.n. 


9-10 





e\e' 


e\e° 


- 


11-12 





e\e' 


e\e\e^e' 


n.n.n. 


13-14 


p3 


£lO £ll £03 £33 


£lO £21 £02 £03 £32 £33 


n.n. 


15-16 


p3 


e\e\e\e' 


e°,e\e\e' 


n.n. 


17-18 





£lO £ll £03 £33 


£lO £ll £02 £32 


n.n. 


19-20 





e\e\e\e' 


e°,e' 


n.n. 


21-22 





£10 £21 £22 £23 


£lO £21 £22 £23 


n.n.n. 


23-24 





e\e\e\e^ 


e',e\e',e' 


n.n. 


25-26 





e\e\e\e' 


e',e',e\e' 


n.n. 


27-28 





£ll £l2 £l3 £20 


£lO £ll £l3 £22 


n.n. 


29-30 





e\e\e\e' 


e\e',e^e' 


n.n. 



pseudo-spin rotation symmetry. However, because of the un- 
usual transformation property of the non-Kramers pseudo- 
spins under time reversal, the operators corresponding to 
Qi? , Qf ? , Jh , Jfj are odd under time reversal, a symmetry of 
the pseudo-spin Hamiltonian. Hence if any of the above op- 
erators gain a non-zero expectation value in the ground state, 
then the corresponding spin-orbital liquid breaks time rever- 
sal symmetry. While this can occur in principle, we check 
explicitly (see Appendix C) that in all the spin-orbital liquids 
discussed above, the expectation values of these operators are 
identically zero. This provides a non-trivial consistency check 
on our PSG calculations. 

We now briefly dicuss the effect of the fluctuations about 
flie mean-field states. In the absence of pairing channels (both 
singlet and triplet) the gauge group is f7(l). In this case, the 
fluctuations of the gauge field about the mean field (Eq. 15) 
are related to the scalar pseudo-spin chirality ^i • 6*2 x 6*3, 
where the three sites form a triangle.'^ Such fluctuations are 
gapless in a C/(l) spin liquid. It is interesting to note that the 
scalar spin-chirality is odd under time-reversal symmetry and 
it has been proposed that such fluctuations can be detected 
in neutron scattering experiments in presence of spin rotation 
symmetry breaking.-" In the present case, however, due to the 
presence of spinon pairing, the gauge group is broken down 
to Z2 and the above gauge fluctuations are rendered gapped 
through Anderson-Higg's mechanism. '^^ 

In addition to the above gauge fluctuations, because of the 



triplet decouplings which break pseudo-spin rotational sym- 
metry, there are bond quadrupolar fluctuations of the pseudo- 
spins Q" (Eq. 35), as well as vector chirality fluctuations Jij 
(Eq. 36)^''^^ on the bonds. These nematic and vector chirality 
fluctuations are gapped because the underlying pseudo-spin 
Hamiltonian (Eq. 4) breaks pseudo-spin-rotation symmetry. 
However, we note that because of the unusual transformation 
of the non-Kramers pseudo-spins under time reversal (only 
the z— component of pseudo-spins being odd under time re- 
versal), Q}? OS^.Jh and J^^ are odd under time reversal. 
Hence, while their mean field expectation values are zero (see 
above), the fluctuations of these quantities can in principle lin- 
early couple to the neutrons in addition to the z— component 
of the pseudo-spins. 



Having identified the possible Z2 spin-orbital liquids, we 
can now study typical dynamic structure factors for these 
spin-orbital liquids. In the next section we examine the typi- 
cal spinon band structure for different spin-orbital liquids ob- 
tained above and find their dynamic spin structure factor. 




M' 



K' 



FIG. 3: The spin structure factor for an ansatz in spin liquid 17, with 
the spin variables transforming as a &amers doublet. 




M' 



K' 



FIG. 4: The spin structure factor for an ansatz in spin liquid 17, with 
the spin variables transforming as a non-Kramers doublet. 



IV. DYNAMIC SPIN STRUCTURE FACTOR 

We compute the dynamic spin structure factor 
dt 



S{q,uj) 



27r 



E' 



,iq-{ri-rj 



^ «(t)a^no)), 

a=l,2,3 

(37) 



for an example ansatz of our spin liquid candidates, in order to 
demonstrate the qualitative differences between the Kramers 
and non-Kramers spin-orbital liquids. In the above equation, 
the pseudo-spin variables are defined in a global basis (with 
the z-axis perpendicular to the Kagome plane). In computing 
the structure factor for the non-Kramers example, we include 
only the a^ components of the pseudo-spin operator in the lo- 
cal basis, since only the z-components carry magnetic dipole 
moment (see discussion before). Hence, only this component 
couples linearly to neutrons in a neutron scattering experi- 
ment. 

Eq. 37 fails to be periodic in the first Brillouin zone of 
the Kagome lattice'*, as the term r.i — tj in eq. 37 is a half- 
integer multiple of the primitive lattice vectors when the sub- 
lattices of sites i and j are not equal. As such, we examine the 
structure factor in the extended brillouin zone, which consists 
of those momenta of length up to double that of those in the 



first brillouin zone. We plot the structure factor along the cut 
r -^ M' -^ K' ~^T, where M' = 2M and K' = 2K. We 
examine the structure factors for two ansatz of spin liquid # 
17 which has both Kramers and non-Kramers analogues. 

As expected, we find that the structure factor has greater 
weight in the case of a Kramers spin liquid. This is par- 
tially due to the fact that the moment of the scattering par- 
ticle couples with all components of the spin, rather than sim- 
ply the z-component. In addition, we note that the presence 
of terms allowed in the non-Kramers spin-orbital liquid in- 
duce the formation of a gap, which is absent for the Kramers 
case with up to second nearest neighbour singlet and triplet 
terms in this particular spin-orbital liquid. Qualitative and 
quantitative differences such as these, which can be observed 
in these structure factors between Kramers and non-Kramers 
spin-orbital liquids, provides one possible distinguishing ex- 
perimental signature of these states. We shall not pursue this 
in detail in the present work. 



V. DISCUSSION AND POSSIBLE EXPERIMENTAL 
SIGNATURE OF NON-KRAMERS SPIN-ORBITAL LIQUIDS 

In this work, we have outlined the possible Z2 spin-orbital 
liquids, with gapped or gapless fermionic spinons, that can be 
obtained in a system of non-Kramers pseudo-spin- l/2s on a 
Kagome lattice of Pr+'^ ions. We find a total of thirty, 10 more 
than in the case of corresponding Kramers system, allowed 
within PSG analysis in presence of time reversal symmetry. 
The larger number of spin-orbital liquids is a result of the dif- 
ference in the action of the time-reversal operator, when real- 
ized projectively. We note that the spin-spin dynamic struc- 
ture factor can bear important signatures of a non-Kramers 
spin-orbital liquid when compared to their Kramers counter- 
parts. Our analysis of the number of invariant PSGs leading to 
possibly different spin-orbital liquids that may be realizable in 
other lattice geometries will form interesting future directions. 

We now briefly discuss an experiment that can play an im- 
portant role in determining non-Kramers spin-orbital liquids. 
Since the non-Kramers doublets are protected by crystalline 
symmetries, lattice strains can linearly couple to the pseudo- 
spins. As we discussed, the transverse {x and y) components 
of the pseudo-spins {a^, a^} carry quadrupolar moments and 
hence are even under the time reversal transformation. Fur- 
ther, they transform under an Eg iiTeducible representation 
of the local D^d crystal field. Hence any lattice strain which 
has this symmetry can linearly couple to the above two trans- 
verse components. It turns out that in the crystal type that 
we are concerned, there is indeed such a mode related to the 
distortion of the oxygen octahedra. Symmetry considerations 
show that the Unear coupUng is of the form Egicr^ + Eg2(J^ 
{{Egi,Eg2} being the two components of the distortion in 
the local basis). The above mode is Raman active. For 
a spin-liquid, we expect that as the temperature is lowered, 
the spinons become more prominent as deconfined quasipar- 
ticles. So the Raman active phonon can efficiently decay into 
the spinons due to the above coupling channel. If the spin 
liquid is gapless, then this will lead to anomalous broaden- 
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ing of the above Raman mode as the temperature is lowered, 
which, if observed, can be an experimental signature of the 
non-Kramers spin-orbital liquid. The above coupling is for- 
bidden in Kramers doublets by time-reversal symmetry and 
hence no such anomalous broadening is expected. 



Acknowledgments 



By requiring that this matrix be invariant under this transfor- 
mation, we can see that this potential only contains matrix 
elements for mixing of states which have the z-component of 
angular momentum which differ by 3. Thus, our eigenstates 
are mixtures of the \m = 4), \m = 1), and \m — —2) states, 
of the \m = 3), \m = 0), and \m = —3) states, and of the 
\m = —4), \m = —1), and \m — 2) states. 
In addition to this, we have the transformation properties 
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and 



r|m) = (-ir|-m) 



a\m) = (-l)™|-m) 



(A7) 



(A8) 



Appendix A: Crystal Field Effects 

In this appendix, we explore the breaking of the J = 4 spin 
degeneracy by the crystalline electric field. The oxygen and 
TM ions form a D^d local symmetry environment around the 
Pj,3+ ions, splitting the ground state degeneracy of the elec- 
trons. This symmetry group contains 6 classes of elements: E, 
2C3, 3C2, i, 25*6, and icjd, where the C3 are rotations by 27r/3 
about the local z axis, the C2 are rotations by vr about axis per- 
pendicular to the local z axis, i is inversion, 5*6 is a rotation 
by 47r/3 combined with inversion and ad is a reflection about 
the plane connecting one corner and the opposing plane, run- 
ning through the Pr molecule about which this is measured 
(or, equivalently, a rotation about the x axis combined with 
inversion). For our J=4 manifold, these have characters given 
by 



r'\E) 



X 



(4), 



2*4 + 1 



sin(37r) 
sin(7r/3) 



x''Had) = x''\.) = '^^^^^^^-^--^'^ 



sin(7r/2) 



= 1 = X 



x'-'Ks,) 

ic'2) 



(Al) 
(A2) 

(A3) 



where the latter equalities are given by the fact that our J=4 
manifold is inversion symmetric. Thus, decomposing this in 
terms of D^d irreps, our 1=4 manifold splits into a sum of 
doublet and singlet manifolds as 



r 



i=4 



3E, + 2Au 



i2ff- 



(A4) 



To examine this further, we need to consider the matrix ele- 
ments of the crystal field potential between the states of differ- 
ent angular momenta. We know that this potential must be in- 
variant under all group operations of D^d, so we can examine 
the transformation properties of individual matrix elements, 
{m\V\m'). Under the C3 operation, these states of fixed m 
transform as 



CsH 



[u 



e 3 



(A5) 



and thus the matrix elements transform as 



C3 : {m\V\m') -^ (m|(C3)- VC3|m') = a;'"'-"(m|y|m'). 

(A6) 



(where the operators for time reversal and reflection are 
bolded for future clarity). Inversion acts trivially on these 
states, as we have total angular momentum even. Thus our 
time-reversal and lattice reflection (about one axis) symme- 
tries give us doublet states of eigenstates a\m — 4) + (3\ni — 
1) — 7|to = —2) and a\m — —4) — /3|m = —1) — 7|m — 2) 
(with a, /3, 7 e 5R in order to respect the time reversal symme- 
try) for the three eigenstates of V in these sectors. The eigen- 
states of the \m — 3), \m = 0), and |to = —3) portion of V 
must therefore spUt into three singlet states, by our represen- 
tation theory argument A4. Due to the expected strong Ising 
term in our potential, we expect the eigenstate with maximal J 
to be the ground state, meaning that to analyze the properties 
of this ground state we are interested in a single doublet state, 
one with large a (close to one). We will restrict ourselves 
to this manifold from this point forward, and define the two 
states in this doublet as 

1+) =a|m = 4) +/3|to= 1) -7|m = -2) (A9) 

|-) = a\m = -4) - f3\m = -1) - 7|m = 2). (AlO) 

We shall also refer to states of angular momentum \m = n) 
as \n) for simplicity of notation. 



Appendix B: Gauge transformations 

We begin by describing the action of time reversal on our 
ansatz. The operation is antiunitary, and must be combined 
with a spin transformation a^ in the case of non-Kramers dou- 
blets. As a result, the operation acts as T : ^" E^F'^ — > 

^," S^I]"*S^r'^*. However, we can simplify this consid- 
erably by performing a gauge transformation in addition to 
the above transformation, which yields the same transforma- 
tion on any physical variables. The gauge transformation we 
perform is iT^, which changes the form of the time reversal 



operation to T : Q^Y.'^VP -^ ^'^P*y,^y.°'*y.^y'^y(^*T'^ = 
I'^^T.^'VP, where S,"^ = ^"^^ if a e {1, 2} and (°'l^ = ~("f^ if 
ae {0,3}. 

On the Kagome lattice, the allowed form of the gauge trans- 
formations has been determined by Yuan-Ming Lu et al. '^ For 
completeness, we will reproduce that calculation, valid also 



for our spin triplet ansatz, here. The relations between the 
gauge transformation matrices. 



[GT{i)f=r,Tl, 

G^^^i^)GU^)GTA^)GTiT^\^)) = vt.tI, 
GUi)GU^)G„{^)GT{<J-Hi)) ^ V.tI, 
Gl^{^)GU^)Gs,i^GT{S^\^)) - t^s^tI, 
4jTf 1W)4^(*)Gt,(0Gt,(T2-1(j)) - ,712/, 
GseiSeH^))GseiSe\^))GsASliz)) 
xGsASmGseiS5{^))Gs,{^) = VsJ, 
Gi{T-\^))G!,A^)GA^)GTA<^{^)) = v.tJ, 

Gl{T^\i))G^^A^G,^)GTA<y{i)) = VaTj, 
GUSe{^))Gs,{Se{^))G,{^)Gs,{<J{^)) = ii„sj, 

GliT,-\z))Gli^)GsA^)GTASe\^)) = Vs^tJ, 
Gl^{T^^T,{^)G^TSTi{i))GTATi(i)) 



Gs,{^)GTASe'i^))^Vs, 



'aT2-^ 1 



(Bl) 
(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 

(BID 

(B12) 

(B13) 



are valid for our case as well, due to the decoupling of spin 
and gauge portions of our ansatz. In the above, the relations 
are valid for all lattice sites i — (a;, y, s), I is the 4x4 identity 
matrix, and the G5 matrices are gauge transformation matri- 
ces generated by exponentiation of the F matrices. The ry's 
are ± 1 , the choice of which characterize different spin liquid 
states. In deriving this form of the commutation relations, we 
have included a gauge transformation iV"^ in our definition of 
the time reversal operator, as this simplifies the effect of the 
operator on the mean field ansatz. 

We turn next to the calculation of the gauge transforma- 
tions. We look first at the gauge transformations associated 
with the translations. We can perform a site dependent gauge 
transformation W{i), under which the gauge transformations 
associated with the translational symmetries transform as 

GTA^)^W{i)GTAi)WHi-x) (B14) 

GTAi)^W{i)GTAi)WAi-y). (B15) 

As such, we can choose a gauge transformation W(i) to sim- 
plify the form of Gt^ and Gt2 ■ Using such a transformation, 
along with condition B7, we can restrict the form of these 
gauge transformations to be 



GrAi) =rfi2l GtA^)^I. 



(B16) 



To preserve this choice, we can now only perform gauge 
transformations which are equivalent on all lattice positions 
{W{x,y,s) — W{s)) or transformations which change the 
shown matrices by an IGG transformation. 

Next, we look at adding the reflection symmetry a. Given 
our formulae for Gti and Gt2, along with the relations be- 



tween the gauge transformations, we have that 

GUT2Ai))G.ii)rif2^VaTj (B17) 

GUTAA^))Ga{^vf2 = VaT^I- (B18) 

Defining Gct(0,0,s) = gcr(s), we have, by repeated application 
of the above. 



Next, using 



G„{x, y, s) = vlT^Vt2VaT29a{s). 

G^{a{i))G^{i)^7^„I 



we find that 



'qal = G„(y, X, a{s))Gcr{x, y, s) 

= {71aT,VaT2f'^'" 9a{cr{s))ga{s). 



(B19) 
(B20) 

(B21) 

(B22) 
(B23) 



Since this is true for all x and y, 7]aTi''la-T2 ~ 1 and thus 
VaTi = JhT2 and ga{cr{s))ga{s) = v^I (where a{u) = 
u, a{v) — w and <t{w) — v). Our final form for the gauge 
transformation is 



Gcrix,y,s) = ?7aT; ^125^^(5) ■ 



(B24) 



Next we look at adding the Sq symmetry to our calcula- 
tion. We can do an IGG transformation, taking GTi{Ti{i)) 
to riS(iT2GTi{Ti{i)), with the net effect being that riSgT2 be- 
comes one (previous calculations are unaffected). We now 
have that 



GlST2'TA^))GsJtH2^I 



(B25) 
{s = u,v) (B26) 



GlAT2\^))GsA^)ri^r' = VSeTj 

GlAT2\i))GsA^)Vi2=Vs,Tj {s = w). (B27) 
Defining Gs^ (0, 0, s) = gs^ (s), we find that 



r~< ( \ n(n~l)/2 i \ 



rt I \ x{x — \')l2+y-\-xy x+y i \ 



'^Sa[.X,y,S} — ri^2 VSeTidSeKS) 



(B28) 

(s — u,v) 
(B29) 

(s = w). 

(B30) 



Using the commutation relation between the a and Sq gauge 
transformations, we find that 

■>l<ysj = vlTM2VhTi9Uv)gse{v)9a{u)gse{u) (B31) 

= vlTM2VhTiVaga{w)gse{v)ga{u)gse{u) (B32) 

giving us that ?7o-Ti'7i2»7S6Ti = 1 and 

gaiu)gseiu)gaiw)gse{v) = VaSeVal- A. similar calcu- 
lation on a different sublattice gives us 

r]asj = VanVl2Vlen9l{w)gSeMgaiv)gSe{w) (B33) 

= Van Vl2'n%T^ Verger {v)gSe {w)ga {v)gSe {w) (B34) 
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TABLE II: We list the solutions of Eq. B43 - B54, along with a set of gauge transformations which realize these solutions. 



No. 


VT Vt 


T VSeT 


Vt Vt 


Se VSe 


Vl2 


9t(u) 


9t{v) 


9t{w) 


gse (m) 


9Se {'") 


9SeM 


1,2 


-I 1 


1 I 1 


±1 


±1 


pO 


pO 


pO 


pO 


pO 


pO 


3,4 


-I 1 


1 I 


1 Tl 


±1 


pO 


pO 


pO 


pO 


-r" 


ir' 


5,6 


-I 1 


-1 


I 


1 Tl 


±1 


pO 


pO 


pO 


ir^ 


iV' 


it' 


7,8 


-I 1 


1 -1 


1 =Fl 


±1 


^V^ 


pO 


-r« 


pO 


^V^ 


pO 


9,10 


-I 1 


1 -1 1 


±1 


±1 


ir^ 


pO 


-r" 


pO 


-iV^ 


ir' 


11,12 


-I 1 


-1 


-1 1 


=Fl 


±1 


iV^ 


pO 


-r" 


iV' 


-iV^ 


ir' 


13,14 


-I 


1 -1 


-1 


1 Tl 


±1 


it' 


ir^ 


iV' 


iV' 


iV' 


it' 


15,16 


-I 


11-11 


±1 


±1 


iV' 


iV' 


iV' 


pO 


pO 


pO 


17,18 


-I 


11-11 


Tl 


±1 


jr» 


it' 


iV' 


pO 


pO 


ir' 


19,20 


-I 


1 -1 


-1 1 


Tl 


±1 


iV' 


iV' 


iV' 


iV' 


-it' 


ir' 


21,22 




1 I 1 


±1 


±1 


pO 


pO 


pO 


pO 


pO 


pO 


23,24 




1 I 


1 Tl 


±1 


pO 


pO 


pO 


pO 


-r" 


ir' 


25,26 




1 -1 


1 Tl 


±1 


ir^ 


pO 


-r" 


pO 


iV' 


pO 


27,28 




1 -1 1 


Tl 


±1 


iV' 


pO 


-r« 


pO 


-it' 


ir' 


29,30 




1 -1 1 


±1 


±1 


it' 


pO 


-r" 


pO 


-it' 


it' 



giving us {ga{v)gse{w))'^ = V-jSeVal- A Z2 (IGG) gauge 
transformation of the form W{x, y, s) — v^t changes 770.^1 
to 1 . Using the cyclic relation of the gauge transformations 
related to the Sq operators, we find 



VSel = Vl2{9Seiw)gSeiv)gSeiu)? 



giving us that 



[9Se{w)gSeiv)9Se{u)? = VSeVl2l- 



(B35) 



(B36) 



Next we turn to the time reversal symmetry. Similar meth- 
ods to the above give us that 



[Grii]? - VtI 



(B37) 
(B38) 
(B39) 



The first of these relations tells us that Grii) is either the 
identity (for vt = 1) or id ■ k (for vt — —1, where \d\ — 1. 
Defining GT(0,0,s)=5r(s), 



Grix^y^s) = VtitVt2t9t{s) 



(B40) 



and further, using the commutation relations between the a 
and T gauge transformations and the 5*6 and T gauge trans- 
formations, 

9Us)gUs)gAs)gTHs))VTtTVTtT = V.tI (B41) 

9leis)9ks)9Se{s)gTiS,\s))4\^^\^^J^^ = VSeTl-iBAl) 



Because this is true for all x and y, and fi{i) is not equal to 
/2(*)' VTiT — 11T2T — 1- If Gt(*) = ia ■ K, we perform a 
gauge transformation W on Grii) such that W^GT{'i)W = 
in^ (as this is the same on all sites, it does not affect our gauge 
fixing for the translation gauge transformations). Collecting 
the necessary results for further use, 



(B43) 
(B44) 
(B45) 

s ^u,v (B46) 

s^w (B47) 



GT2{x,y,s) =/ 
Ga{x,y,s) = Vi29'7{s) 

Gt{s) = / = gris) m = 1 (B48) 

Gt{s) = iV^ ^ gris) vt = -I (B49) 

g^{ais))g^is) ^ Val (B50) 

ga{v)gse{'u)ga(w)gse{v) = {g„{v)gse{w)f = VaSeVal 

(B51) 

{gSe {w)gSe i.v)gse (w))^ = f?S6'7l2^ (B52) 

g<y{s)gT{(j{s)) = VaT9T{s)g„{s) (B53) 

9Sei.s)gT{S^^{s)) ^vseTgT{s)gse{s)- (B54) 

We also have the gauge freedom left to perform a gauge ro- 
tation arbitrarily at all positions for 77^ = 1 or an arbitrary 
gauge rotation about the x axis for tjt — — 1- 

The solution to the above equations is derived in detail by 
Lu et al}^ and as such we simply list the results in table 11. 
The basic method of obtaining these solutions is as follows: 
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for each choice of Z2 parameter set, we determine whether 
there is a choice of gauge matrices { gs } which satisfy the 
equations B43 - B54. In order to do so, we determine the 
allowed forms of the gs matrices from the equations, then use 
the gauge freedom on each site to fix the form of these. Of 
particular not is the fact that in the consistency equations for 
the g matrices, the terms r\Yi and 775^ only appear multiplied 
together, meaning that for any choice of the gauge matrices 
gs we can choose 7712 = ±1, which fixes the form of 7755. 



where our definition of r\ij differs by a factor of (-1) from that 
of the cited work. We rewrite this in terms of our variables, 
finding 



Appendix C: Relation among the mean-field paramters 



The relation among the different singlet and triplet param- 
eters in terms of ^^ is given by 



4 



Dl = -e 



10 

11 



S27 ' S^ 



4 



*e 



i-'j 

13. 

12 



S^7 I S^ 






m + ^i 



02. 



23 

22. n3 



^3 _ e30 



e3u I e33 






32 



(CI) 



-^7 






(5^'' 



Epiiapva 



E((e 



1:°)^ 



xc-n') 



/tA _»('tOO<jAO _ Y^ <jOa^Aa\ 



(C3) 



Using these, we can derive the form of the bond nematic 
order parameter and vector chirality order parameters, which 
are given in terms of the mean field parameters-' as 

Q^r--l{E^E:;^~h^^'^\E,,{') + h.c. 



D';^D*;'--S^''^\D,,{')+h.c. 



'-'ij ~ n\Xij-^ij Xij ij) 



(C2) 



In particular, we find that J^, J^, Q}^ and Q^^ must be zero 
for all non-Kramers spin liquids, as the terms allowed by sym- 
metry in Eq. 32 and 33 do not allow non-zero values for these 
order parameters. 
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